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Dark energy cosmology is considered in a modified Gauss-Bonnet (GB) model of gravity where 
an arbitrary function of the GB invariant, f(G), is added to the General Relativity action. We show 
that a theory of this kind is endowed with a quite rich cosmological structure: it may naturally lead 
to an effective cosmological constant, quintessence or phantom cosmic acceleration, with a possibility 
for the transition from deceleration to acceleration. It is demonstrated in the paper that this theory 
is perfectly viable, since it is compliant with the Solar System constraints. Specific properties of 
f{G) gravity in a de Sitter universe, such as dS and SdS solutions, their entropy and its explicit 
one-loop quantization are studied. The issue of a possible solution of the hierarchy problem in 
modified gravities is addressed too. 



PACS numbers: 98.70.Vc 

I. INTRODUCTION 



Recent observational data indicate that our uni- 
verse is accelerating. This acceleration is explained 
in terms of the so-called dark energy (DE) which 
could result from a cosmological constant, an ideal 
fluid with a (complicated) equation of state and 
negative pressure, the manifestation of vacuum ef- 
fects, a scalar (or more sophisticated) field, with 
quintessence-like or phantom-like behavior, etc. (For 
a very complete review of a dynamical DE see and 
references therein, for an earlier review, see 0.) The 
choice of possibilities reflects the undisputable fact 
that the true nature and origin of the dark energy 
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has not been convincingly explained yet. It is not 
even clear what type of DE (cosmological constant, 
quintessence or phantom) occurs in the present, late 
universe. 

A quite appealing possibility for the gravitational 
origin of the DE is the modification of General Rel- 
ativity. Actually, there is no compelling reason why 
standard GR should be trusted at large cosmologi- 
cal scales. For a rather minimal modification, one 
assumes that the gravitational action may contain 
some additional terms which start to grow slowly 
with decreasing curvature (of type \/R 0, 0| , In i? 
, Tr 1/i? , string- inspired dilaton gravities , 
etc.), and which could be responsible for the current 
accelerated expansion. In fact, there are stringent 
constraints on these apparently harmless modifica- 
tions of General Relativity coming from precise Solar 
System tests, and thus not many of these modified 
gravities may be viable in the end. In such situa- 
tion, a quite natural explanation for both the cosmic 
speed-up issue and also of the first and second coin- 
cidence problems (for recent discussion of the same. 
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see |g) could be to say that all of them are caused, in 
fact, by the universe expansion itself! Nevertheless, 
one should not forget that some duality exists be- 
tween the ideal fluid equation of state (EoS) descrip- 
tion, the scalar-tensor theories and modified gravity 
|3- Such duality leads to the same FRW dynamics, 
starting from three physically different — but math- 
ematically equivalent — theories. Moreover, even 
for modified gravity, different actions may lead to 
the same FRW dvnamicsfioj. Hence, additional evi- 
dence in favor of one or another DE model (with the 
same FRW scale factor) should be clearly exhibited 
0. 

As a simple example, let us now see how different 
types of DE may actually show up in different ways 
at large distances. It is well known that cold dark 
matter is localized near galaxy clusters but, quite 
on the contrary, dark energy distributes uniformly in 
the universe. The reason for that could be explained 
by a difference in the EoS parameter w — p/p. As 
we will see in the following, the effect of gravity on 
the cosmological fluid depends on w and even when 
— 1 < w < gravity can act sometimes as a repulsive 
force. 

To see the w— dependence on the fluid distribution 
in a quite simple example, we consider cosmology in 
anti-de Sitter (AdS) space, whose metric is given by 



ds^ = dy + e 



2y/l 



Y^{d. 



(1) 



the conservation law of the energy momentum ten- 
sor, V^r^'' gives, by putting v = y, 



dp 1 , 



(2) 



if we assume matter to depend only on the coordi- 
nate y. When w is a constant, we can solve Eq. Q 
explicitly 



po exp 



7^ + ^'^ 



(3) 



Here po is a constant. We should note that > 
when w > or w < —1, and I + 1/w < when 

— l<w<0. Then, for usual matter with w > 0, 
the density p becomes large when y is negative and 
large. In particular, for dust w ~ but w > 0, and 
a collapse would occur. On the other hand, when 

— 1 < w < 0, like for quintessence, the density p 
becomes large when y is positive. In the phantom 
case, with w < —1, the density p becomes large when 



y is negative, although a collapse does not occur. 
When w = —1^ p becomes constant and uniform. 
We may also consider a Schwarzschild like metric: 



ds' 



(4) 



Here (iri^2) expresses the metric of a two-dimensional 
sphere of unit radius. Then, the conservation law 
VpT'"' with V = r gives, 

dp dv , sr. 

-T + ^(P + P) = ^ ■ (5) 
dr dr 

When w is constant, we can solve (O and obtain 



Po exp - IH ]^ ■ (6) 

















w J 



Here pq is a constant again. In particular, in the 
case of the Schwarzschild metric, 



with horizon radius tq, we flnd 



P = Po 



1 



r 



(7) 



(8) 



Then, when w > oi w < — 1, pisa decreas- 
ing function of r, that is, the fluid is localized near 
the horizon. Specifically, in the case of dust with 
w ~ Q, the fiuid collapses. On the other hand, when 
— I<ii'<0,pisan increasing function of r, which 
means that the fluid delocalizes. When w = 0, the 
distribution of the fluid is uniform. 

The above results tell us that the effect of grav- 
ity on matter with — 1 < w < is opposite to 
that on usual matter. Usual matter becomes dense 
near a star but matter with — 1 < w < becomes 
less dense when approaching a star. As is known, 
cold dark matter localizes near galaxy clusters but 
dark energy distributes uniformly within the uni- 
verse, which would be indeed consistent, since the 
EoS parameter of dark energy is almost — 1 . If dark 
energy is of phantom nature {w < —1), its density 
becomes large near the cluster but if dark energy is 
of quintessence type {—l<w< —1/3), its density 
becomes smaller. 

In the present paper the (mainly late-time) cos- 
mology coming from modified Gauss-Bonnet (GB) 
gravity, introduced in Ref.|llj|. is investigated in de- 
tail. In the next section, general FRW equations 
of motion in modified GB gravity with matter are 
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derived. Late-time solutions thereof, for various 
choices of the function f{G), are found. It is shown 
that modified GB gravity may indeed play the role 
of a gravitational alternative for DE. In particu- 
lar, we demonstrate that this model may naturally 
lead to a plausible, effective cosmological constant, 
quintessence or a phantom era. In addition, /(G) 
gravity has the possibility to describe the inflation- 
ary era (unifying then inflation with late-time accel- 
eration), and to yield a transition from deceleration 
to acceleration, as well as a natural crossing of the 
phantom divide. It also passes the stringent solar 
system tests, as it shows no correction to Newton's 
law in flat space for an arbitrary choice of /(G), 
as well as no instabilities. Sect. 3 is devoted to 
the study of the de Sitter universe solution in such 
model. The entropies of a SdS black hole and of a dS 
universe are derived, and possible applications to the 
calculation of the nucleation rate are discussed. In 
Sect. 4, the quantization program at one-loop order 
for modified GB gravity is presented. This issue is 
of the essence for the phantom era, where quantum 
gravity effects eventually become important near the 
Big Rip singularity. Sect. 5 is devoted to the gener- 
alization of modified gravity where F = F{G,R). 
This family of models looks less attractive, given 
that only some of its specific realizations may pass 
the solar system tests. Nevertheless, it can serve to 
discuss the origin of the cosmic speed-up as well as 
a possible transition from deceleration to accelera- 
tion. In Sect. 6, the important hierarchy problem 
of particle physics is addressed in the framework of 
those modified gravity theories. It is demonstrated 
there that this issue may have a natural solution in 
the frame of J-{R) or F{G, i?)-gravity. The last sec- 
tion is devoted to a summary and an outlook. In 
the Appendix, an attempt is made to construct zero 
curvature black hole solutions in the theory under 
discussion. 



II. LATE-TIME COSMOLOGY IN 
MODIFIED GAUSS-BONNET GRAVITY 

Let us start from the following, quite general ac- 
tion for modified gravity [isj 

S = J d^x^g (/ {R, Rf^^R'^", R^.p^R'^"'"') + L„ 

(9) 

Here £,„ is the matter Lagrangian density. It is 
not easy to construct a viable theory directly from 
this general class, which allows for non-linear forms 



for the action. One must soon make use of sym- 
metry considerations, which lead to theories which 
are more friendly, e.g., to the common Solar System 
tests. Specifically, we shall restrict the action to the 
following form: 



S = d":. 



(10) 



(11) 



-g{F{G,R)+Cra) ■ 
Here G is the GB invariant: 

Varying over g^^: 

O^T'"' + ^g>"'F{G,R) 

-2Fg{G, R)RR'"' + AFg{G, R)R''pR''P 
-2Fg{G, R)R^p''^R\,, - 4^^g(G, R)R'''"'''Rp^ 
+2 {W^WFaiG, R)) R - 2g>"' {W'^Fg{G, R)) R 
-4 (VpV''Fg(G, R)) R"" - 4 {WpW^FaiG, R)) R^p 
+4: {V^Fg{G, R)) R^'' + ig'"' (VpV,Fg(G, R)) RP" 
-4(VpV„FG(G,i?))i?'^'"''^ 
-Fr{G, R)R^"' + V^'VFRiG, R) 
-gP^V^FniG^R) , (12) 

TP'^ being the matter-energy momentum tensor, and 
where the following expressions are used: 



dG 



dR 



(13) 

The spatially-flat FRW universe metric is chosen as 



ds" = -de + a{tfY,{dxy 



(14) 



i=l 



Then the (t, t)— component of H12|) has the following 
form: 



= GFg{G,R)- F{G,R)-2AH' 
'dH 



^dFG{G,R) 



dt 

dFniG, R) 



+ P71 



(15) 



where pm is the energy density corresponding to 
matter. Here, G and R have the following form: 



G = 2a(^H^H + H^^ , R = q[h + 2H^^ . (16) 
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In absence of matter = 0), there can be a de 
Sitter solution (H = Hq — constant) for (fTK|l . in 
general (see 0). One finds Hq by solving the alge- 
braic equation 

= 24H^Fg{G, R) - F{G, R) + 6H^Fr{G, R) . 

(17) 

For a large number of choices of the function 
F(G, R), Eq. 1201) has a non-trivial {Hq ^ 0) real so- 
lution for Hq (the de Sitter universe) . The late-time 
cosmology for the above theory without matter has 
been discussed for a number of examples in refs.|llj. 
In this section, we restrict the form of F{G, R) to 

be 



FiG, R) 



1 
2k 



,R + fiG) , 



(18) 



where = SttGn, Gm being the Newton constant. 
As will be shown, such an action may pass the So- 
lar System tests quite easily. Let us consider now 
several different forms of such action. By introduc- 
ing two auxilliary fields, A and B, one can rewrite 
action itTUI) with l(TH)l as 



S = 



j d^x^giJ-^R + B{G-A) 



(19) 



Varying over i?, it follows that A ^ G. Using this 
in ((T^ . the action (^Ul with ((TH|l is recovered. On 
the other hand, varying over A in (|19|) . one gets 
B = f'{A), and hence 

1 



S = / d^x^i^^R + f'iA)G 

-Af'{A) + fiA)). (20) 



By varying over A, the relation A = G is obtained 
again. The scalar is not dynamical and it has no 
kinetic term. We may add, however, a kinetic term 
to the action by hand 



S 



-d^Ad^A 



+f{A)G - Af'{A) + f{A) 



(21) 



Here e is a positive constant parameter. Then, one 
obtains a dynamical scalar theory coupled with the 
Gauss-Bonnet invariant and with a potential. It is 
known that a theory of this kind has no ghosts and 
it is stable, in general. Actually, it is related with 
string-inspired dilaton gravity, proposed as an alter- 
native for dark energy |2i]. Thus, in the case that 



the limit e ^ can be obtained smoothly, the cor- 
responding /(G) theory would not have a ghost and 
could actually be stable. This question deserves fur- 
ther investigation. 

We now consider the case ^ 0. Assuming that 
the EoS parameter w = Pm/pm for matter (pm is 
the pressure of matter) is a constant, then, by using 
the conservation of energy: pm -I- 37? {pm + Pm) = 0, 



we find p — 
is given by 



poa 



-3{l+w} 



We also assume that /(G) 



/(G) = /o \Gf 



(22) 



with constants /o and /3. If /3 < 1/2, the /(G) term 
becomes dominant, as compared with the Einstein 
term, when the curvature is small. If we neglect 
the contribution from the Einstein term in (|15|l with 
H18|) . assuming that 



ao (ts - t) 



when Hq > 
when ho < 



(23) 



the following solution is found 



ho = 



ao 



4/3 



- /o(/3"l) 

- (/lo-l)/90 

x{ho-l + 4/3) 



{24\hl i^l + ho)\Y 



"3(1 + ™) 



(24) 



One can define the effective EoS parameter Weff as 



WcS 



2H 
3IP 



(25) 



which is less than — 1 if /3 < 0, and for > — 1 as 



Wcff 



2 
3/io 



1 



2/3 



(26) 



which is again less than —1 for /3 < 0. Thus, if 
/3 < 0, we obtain an effective phantom with nega- 
tive ho even in the case when w > —1. In the phan- 
tom phase, a singularity of Big Rip type a.t t — tg 
seems to appear (for the classification of these 
singularities, see 0|). Near this sort of Big Rip sin- 
gularity, however, the curvature becomes dominant 
and then the Einstein term dominates, so that the 
/(G)-term can be neglected. Therefore, the universe 
behaves as a = aot^/'^*-'"''"^-' and, as a consequence, 
the Big Rip singularity will not eventually appear. 
The phantom era is transient. 

A similar model has been found in llq by using a 
consistent version of 1/i^-gravityQ- In general. 
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in the case of jF(i?)-gravity instabilities appear [13 ■ 
These instabihties do not show up for the case of 
/(G)-gravity. 

Note that under the assumption (|23|l . the GB in- 
variant G and the scalar curvature R behave as 

^ _ 2Ahl {ho - 1) 24hl [ho - 1) 

(_T — ;; , or 



t4 



^ ^ 6feo(2feo-l) 6ho (2ho-l) ^^^^ 



t2 



its - ty 



As a consequence, when the scalar curvature R be- 
comes small, that is, when t ov tg — t becomes large, 
the GB invariant G becomes small more rapidly than 
R. And when R becomes large, that is, if t or ts — t 
becomes small, then G becomes large more rapidly 
than R. Thus, if /(G) is given by ^ with /3 < 1/2, 
the /(G)-term in the action H1(J|I with (|18|l becomes 
more dominant for small curvature than the Einstein 
term, but becomes less dominant in the case of large 
curvature. Therefore, Eq. (|24|l follows when the cur- 
vature is small. There are, however, some exceptions 
to this. As is clear from the expressions in H27|) . 
when h{) = —1/2, which corresponds to Woff = —7/3, 
R vanishes, and when /iq = —1, corresponding to 
WeS = —5/3, G vanishes. In both these cases, only 
one of the Einstein and /(G) terms survives. 

In the case when /3 < 0, if the curvature is large, 
the Einstein term in the action (|10|) with (|18|l dom- 
inates and we have a non-phantom universe, but 
when the curvature is small, the /(G)-term domi- 
nates and we obtain an effective phantom one. Since 
the universe starts with large curvature, and the cur- 
vature becomes gradually smaller, the transition be- 
tween the non-phantom and phantom cases can nat- 
urally occur in the present model. 

The case when < /3 < 1/2 may be also con- 
sidered. As l3 is positive, the universe does not 
reach here the phantom phase. When the curva- 
ture is strong, the /(G)-term in the action H1U|) with 
(|18|l can be neglected and we can work with Ein- 
stein's gravity. Then, if w is positive, the matter 
energy density pm should behave as ~ i^^, but 
/(G) goes as /(G) - t""^. Then, for late times 
(large t), the /(G)-term may become dominant as 
compared with the matter one. If we neglect the 
contribution from matter, Eq. (|15|l with (|18|l has 
a de Sitter universe solution where H, and there- 
fore G, are constant, li H — Hq with constant Hq, 
Eq. with lO looks as with lO. As a 
consequence, even if we start from the deceleration 
phase with w > —1/3, we may also reach an asymp- 
totically de Sitter universe, which is an accelerated 



universe. Correspondingly, also here there could be 
a transition from acceleration to deceleration of the 
universe. 

Now, we consider the case when the contributions 
coming from the Einstein and matter terms can be 
neglected. Then, Eq. ((T5|) with l(TS|) reduces to 



= G/'(G) - /(G) - 2AGf"{G)H^ 



(28) 



If /(G) behaves as (1^ . from assumption (|23J), we 
obtain 

0= {(3-1) hi [ho - 1) (/lo - 1 + 4/3) . (29) 
As /iQ = 1 implies G = 0, one may choose 

/lo = 1 - 4/3 , (30) 

and Eq. gives 

2 



WcS 



-1 



3(1 - 4/3) 



(31) 



Therefore, if /3 > 0, the universe is accelerating 
{wcS < —1/3) and if (3 > 1/4, the universe is in 
a phantom phase (woff < ~1)- Thus, we are led to 
consider the following model: 



/(G) - /, \Gf + fi \Gf^ 
where we assume that 



(32) 



(33) 



Here, when the curvature is large, as in the primor- 
dial universe, the first term dominates, compared 
with the second one and the Einstein term, and gives 



-1 > Wcff = -1 



3(1 - 4A) 



> -5/3 . (34) 



On the other hand, when the curvature is small, as 
is the case in the present universe, the second term 
in H32|l dominates, compared with the first one and 
the Einstein term, and yields 



UJeS = -1 



3(1 - 4A) 



< -5/3 . (35) 



Therefore, the theory (|32|l can in fact produce a 
model which is able to describe both inflation and 
the late-time acceleration of our universe in a unified 
way. 



Instead of (|33|l . one may also choose /?/ as 



^ > A > 



(36) 
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which gives 

> Woff > -1 . (37) 

Then, what we obtain is effective quintessence. 
Moreover, by properly adjusting the couplings fi 
and // in (|32|l . we can obtain a period where the 
Einstein term dominates and the universe is in a de- 
celeration phase. After that, there would come a 
transition from deceleration to acceleration, where 
the GB term becomes the dominant one. 

One can consider the system (|32|) coupled with 
matter as in (|15() with H18() . To this end we just 
choose 

^.>l>Pi. (38) 

Then, when the curvature is large, as in the pri- 
mordial universe, the first term dominates, as com- 
pared with the second one and the Einstein term. 
And when the curvature is small, as in the present 
universe, the second term in (|32|l is dominant as 
compared with the first and Einstein's. Then, an 
effective Wcs can be obtained from (|2()|l . In the 
primordial universe, matter could be radiation with 
w = 1/3, and hence the effective w is given by 

2 

w»,cff = -1 + ^ ' (39) 

which can be less than —1/3, that is, the universe 
is accelerating, when /3i > I. On the other hand, 
in the late-time universe matter could be dust with 
w = 0, and then we would obtain 

wi,oS = -l + :^, (40) 
2p/ 

which is larger than 0, if < /?/ < 1/2, or less than 
— 1, if /?; is negative. Thus, acceleration could occur 
in both the primordial and late-time universes, if 

A > 1 , A < . (41) 

Similarly, one can consider DE cosmology for other 
choices of f(G), for instance, InG or other function 
/ increasing with the decrease of G (late universe) . 

Let us address the issue of the correction to New- 
ton's law. Let (7(0) be a solution of (|12|l with (|18|l 
and represent the perturbation of the metric as 
dfj.i' — 5(0)^11' + First, we consider the perturba- 
tion around the de Sitter background which is a so- 
lution of (flTI with fH^ . We write the de Sitter space 



metric as g(o)^tiy, which gives the following Riemann 
tensor: 

R{o)tiupa = H'o {g(a)iip9{o)iya - 9{Q)p(7g(o)pp) ■ (42) 

The flat background corresponds to the limit of 
Hq — > 0. For simplicity, the following gauge con- 
dition is chosen: g^^-^h^^ — V^^^/i^^ = 0. Then 

Eq. (O with (HHI gives 

- ^ {V^h^u - 2H^h^,) + T^, . (43) 

The GB term contribution does not appear except 
in the length parameter l/i?o of the de Sitter space, 
which is determined with account to the GB term. 
This may occur due to the special structure of the 
GB invariant. Eq. (|43|l tells us that there is no cor- 
rection to Newton's law in de Sitter and even in the 
flat background corresponding to — > 0, what- 
ever is the form of / (at least, with the above gauge 
condition). (Note that a study of the Newtonian 
limit in 1/i? gravity (where significant corrections 
to Newton's law may appear), and its extension has 
been done in 0,0|.) For most 1/i? models the cor- 
rections to Newton's law do not comply with solar 
system tests. 

Expression H43|l can be actually valid in the de 
Sitter background only. In a more general FRW uni- 
verse, there can appear corrections coming from the 
/(G) term. We should also note that in deriving 
(|43|) . a gauge condition g^^^-^h^i, = was used, but 
if the mode corresponding to g^Q-^h^^, is included, 
there might appear corrections coming from the 
/(G) term. The mode corresponding to g'^^-^h^i, gives 
an infinitesimal scale transformation of the metric. 
Then, it is convenient to write the metric as 

gpu = e'^<7(o)^^ . (44) 

Here g(o)pD expresses the metric of de Sitter space in 
(|42|) . The Gauss-Bonnet invariant G is correspond- 
ingly given by 

+v2o)9^0a'^0 + 2V(o)^V(o),0a[;)0a('o)0} 

= e-2^ {24i74 -f Vfo) (-12H2V(o)^0 
+29^0V^o)0 - 2V(o)^V(o),(/.V^o)</' 
-l-a^09(^o)<^^(o).0)} ■ (45) 
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The covariant derivative associated with g(o)fii' is 5(0)^"^ 
written here as y{o)fi- The expansion of f{G), with 
respect to (j), is 



•-pa 



= y^{/(24i/o')+2/ (24ij4; 
-48/' {24H^) H^4>^ 

+0 (0^) + total derivative terms} . (46) 



Since the last term contains (^fo)4'^ ? in general, 
there could be an instability. A way to avoid the 
problem is to fine-tune /(G) so that /" [2AH^) van- 
ishes for the solution Hq in (|17|) with p8l) . 



In order to consider a more general case, one ex- 
pands /(G) in the action ({117)1 with lfTH|l as 



/(G) - /(G(o)) + -/"(G(o))<5G2 

+/'(G(0))(52G + O(/i3) , (47) 



6'G = 



+AR''-RP/^h^p 

-4V,V^/i,^ + RP^'P^VpV^K^) 
[h<^n {-R (VV^/i^'^ + V^V/i/) 
-^R'^'^VPVKp -f AR'^'V'^V^h} 



(48) 



(Tfl 



--rc-p- (VpV./i;' - VpV'/i,,,) 

-VV^/iJ^-hV^V^V) 

-Vxh^^V.h"^ + {2h^PV''Kx 
-h^^Vxh^" + h^^VPKx)) 
-AR-- {{V^'Kx + ^uh^ 
-VxhP)V^Kp 

+ (V^/i^;, + V^Kx - Va/i'^J V,/i^^ 
+ (2V^V-Va/i„")V,/i/ 

{2h^PV,Kx 
+ (V^/l^A + V^/l^A - Va V) /i^^)) } 
+i?^'''"^((V,/i,A + V,V 
-Va/i^p) V^V 

+ i'^p.hpx + Vphf^x - '^xhpp) '^uhax 

+Vp((V^/l^A + V^/l^A 



-VA/i..)/i'J}] 



(49) 



Note that the term proportional to SG does not ap- 
pear in (|47|l since the background metric is a solution 
of (|42l) . Here and in the following the index (0) is al- 
ways suppressed, when there cannot be confusion. If 
where G(o) is the Gauss-Bonnet invariant given by we choose the gauge condition V/i^jy = 0, Eqs. H48|l 
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and I49|l have the following form: 

+R.'^''''-VpV,K^) , (50) 
= [hi^n{-R{hP^Rp"' ~h^pRP'^^''^) 

+2R^'' {hp^RP'^ - KpRPD 
+AR-^ {hP^Rp, - hppRP^^P^) 
+2i?';'^"^ {hp^RP, - hppRP/,) 

-RP-'P^WpW.h^} 

+ 1 {-4 ( V^". - h^pR'f.) KRP" 
-h^pRP''^"' - V^h''" - V^V^h^") 
+VpV^Kf, {WPWh^p - VPVPh"" 

-vvhpp + yyph'"')} 

^^{R {2VPKxV^h''p ~ Vxh^ 

-AR'"' {{VPKx + V,hP^ 
-VxhP)V^Kp 

+ [V^hP^ + V^Kx - VxhP^) \I,h\ 

+ (V^ft^A + ^^.Kx - Va V) 
+RP-P'' {{VpKx + y^hpx 
-yxKp)V^Kp 

+ (Vp/ipA + Vp/lpA - ^xhpp) Vyhcx 
+Wp{{W^h„x+^ah^x 

-VxKa)h\)}] . (51) 

Now, we consider the case that H ^ in the FRW 
universe (|14|l . Then, by specifying the dimension, 
the following structure is found 

5G - H'^h + H^V^h , 

+H^hVh+{V^hf . (52) 

For the qualitative arguments that follow, we have 
abbreviated the vector indices and coefficients. Since 
Eq. 1071) contains SG"^ and 5'^G terms, the i/^V^/i- 

term in 5G and (V^/i)^-term in 5^G have a possibil- 
ity to generate the instability. Explicit calculations 



in the FRW universe tell us that the (V^/i) term 
in 5^G vanishes identically, while the H^V^h term 
in 5G has the following form: 

H^y^h term in 5G = -AHV^hu ■ (53) 

For simplicity, we have chosen again the gauge con- 
dition g'l^^-^hpi, = 'VP'hpi, = 0. Then, except for the 

H = Q case, which describes the de Sitter universe, 
there might be an instability. 

Since the JG^-term has a factor /"(G(o)), if one 
properly chooses the form of /(G) and fine-tunse the 
coefficients, it could occur that /" (G(o)) = in the 

present universe. Correspondingly, the term (V^/i)^ 
does not appear in the action and no instability ap- 
pears. 

To summarize, in both cases ()46|l and H52|l . if we 
choose /" = in the present universe, the instability 
does not appear. As an example, one can consider 
the model (jS21)- As 

/"(G) = /.A (A - 1) \Gf'~^+m (A - 1) \Gf'~' , 

(54) 

if we choose the parameters fi, f3i, fi, and j3i to 
satisfy 

= /.A (A - 1) \Gi^)f~^' + fiPi (A - 1) , (55) 

we find /' (G(o)) = and thus there will be no insta- 
bility. In H55() , G(o) is the value of the Gauss-Bonnet 
invariant given by the curvature in the present uni- 
verse. 

One now rewrites (|12|) with p8|l as an FRW equa- 
tion: 

= -^H^ + PG + PM , (56) 
{) = ^{2H + 3i/2) + pG + p,„ . (57) 

Here pa and pQ express the contribution from the 
/(G) term in the action lfTU|l with ((TH|l : 

PG - G/'(G) - /(G) - 24G/"(G)ff3 , 
PG = -Gf{G) + f{G) + 2AGf\G)H^ 
+8G^f"'{G)H^ 

-192/"(G) (^-SH^HH - 6H^H^ - H^H 
-dH^ii - ISH^H^ + AH^H^ . (58) 

One can view the contribution from the /(G) term as 
a sort of matter satisfying a special (inhomogeneous) 
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EoS 13 iof usual EoS with time-dependent bulk 
viscosity l20j) of the form 

- pG+PG~8G^f"'{G)H^ 

+192/"(G) ( 



(59) 

In particular, in the case of de Sitter space, owing 
to the fact that the Hubble rate H, and therefore G, 
are constant, we find = + Pg- In the case that 
f{G) is given by (|22|l and one further assumes H23|) . 
one gets 

/lo - 1 + 4/3 



PG 

PG 



fo\GfW-l) 



hn-l 



fo\Gf{f3-l) 

3feg - (3 + 8l3)ho + 16(3^ - 4/3 
^ 3ho {ho - 1) ■ ^^^^ 

It follows that the effective EoS wq = Pg/pg for the 
/(G) part is given by 

3/i2 » (3 + 8/3)/io + 16/32 - 4/3 

= , , • (61) 



In absence of matter (p. 
may be rewritten as 



3ho {ho-l + 4/3) 

Pra 



0), Eq. (EH 



^ -^PG , H = — -(pg+Pg) ■ 
b 2 

Then, by using the expression H16|l . we find 
G = — —Pg C^Pg + 3pg) , 



(62) 



(63) 



k'^P 



G = i^K'Hp + K,''ipG+PG){pG + 5pG)]J-f ^ 

and by using the first equation in (|58|1 . the effective 
equation of state is: 

= -pG 

2^4 / 2^4 
^PG C^PG + 3pg) f ( ^PG C^PG + Spg) 

f 2^4 

-/ ( — ^PG (2pG + 3pg) 

5/2 



-(24)2<i//(:^PG 



^2 /^2 

(PG + Pg) (pg + 3pg) ( yPG 



X/" 



2^4 



PG (2pG + 3pg) 



(64) 



which has the form oi F {^pg.,Pg,H^ ==0 1^. At 
the dynamical level, this demonstrates the equiva- 
lency between modified GB gravity and the effective 
inhomogeneous EoS description. It may be of inter- 
est to study the choice of /(G) which leads to the 
inhomogeneous generalization of EoS p ~ —p — Ap°' 
suggested and studied in |^ as this may be easily 
compared with standard ACDM cosmology. 

Summing up, late-time cosmology in modified GB 
gravity with matter was studied. It has been shown 
that effective DE of quintessence, phantom or cos- 
mological constant type can be actually produced 
within such theory, with the possibility to unify it 
also with primordial inflation. Moreover, the tran- 
sition from the deceleration to the acceleration era 
easily occurs in some versions of the /(G) theory 
which comply with the solar system tests (no in- 
stabilities appear, no corrections to Newton's law 
follow) . 



III. DE SITTER SOLUTION IN MODIFIED 
GAUSS-BONNET GRAVITY MODELS 

In this section, we will investigate the properties 
of some black hole solutions within the modified GB 
gravity scheme. The role played by static, spheri- 
cally symmetric solutions, as the Schwarzschild one, 
with regard to solar system tests, is well known. Due 
to the absence of this kind of zero-curvature BH in 
the /(G) theory (see the Appendix) deSitter space 
as well as SdS BH need her special consideration, 
what it is done below. 

For the general model, it can be shown that the 
equations of motion have the following structure (no- 
matter case) 
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F(G, R)g^iy — Fii{G, R)Rfj,i, 



pap 



— '2'Fg{G, R) [RR^u — "^RfipR^ + R/ipapRJ" 

- 2R^^,pR'"') +E^,{V^Fg{G,R),VMG,R)) . 

where S^j^ is vanishing when G and R are constant. 

In this section we are interested in finding con- 
dition assuring the existence of solutions of the de 
Sitter type (including SdS BH). In such case, the 
Gauss-Bonnet invariant and the Ricci scalar are con- 
stant and the Gauss-Bonnet invariant reads 



R — i?o G — Go — 77^0 ■ 



(65) 



Assuming the maximally symmetric metric solution. 
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one gets 

'^Fr{G,R)R^^^^, = {F{Gq,Ro) — GoFcoiGo, Ro)) g1i,i£i6) 

Note that if F{G, R) = R + aG, namely there is 
only a linear term in G, one gets the ordinary Ein- 
stein equation in vacuum, as it should be, because 
in this case, F{G, R) contains the Hilbert-Einstein 
term plus a topological invariant. 

Taking the trace of H66(l . the condition follows 



RoPb.o (Go, i?o) 



[F{Gq,Ro) — GoFg„{Go, Ro)]{B7) 



This condition will play an important role in the 
following and it is equivalent to the condition H17|) . 
As an example, when F{G,R) = R + f{G), (here 
= 1), one has 



Go/' (Go) -/(Go) = 



(68) 



In general, solving Eq. (|67|l in terms of Rq, one can 
rewrite the maximally symmetric solution as 



pO _ _ A „o 



(69) 



which defines an effective cosmological constant. For 
example, when /(G, R) = R + f{G), one has 



Ae// = I (Go/'(Go) - /(Go)) 

Thus, if /(G) = -aG'^, one has, 

P 



2a(l-/3)(i 



— Hq 



When (3 is small, a > 0, one obtains Rq ^ 
with the choice /3 = —1/2, a > 0, one has 

Rn = 6i/*V3a, 



(70) 

(71) 
a, while 

(72) 



and the corresponding effective cosmological con- 
stant reads 



(73) 



As in the pure Einstein case, one is confronted 
with the black hole nucleation problem 1221. We 
review here the discussion reported in refs. |22ll2^ . 

To begin with, we recall that we shall deal with 
a tunneling process in quantum gravity. On general 
backgrounds, this process is mediated by the asso- 
ciated gravitational instantons, namely stationary 



solutions of Euclidean gravitational action, which 
dominate the path integral of Euclidean quantum 
gravity. It is a well known fact that as soon as 
an imaginary part appears in the one-loop parti- 
tion function, one has a metastable thermal state 
and thus a non vanishing decay rate. Typically, this 
imaginary part comes from the existence of a neg- 
ative mode in the one-loop functional determinant. 
Here, the semiclassical and one-loop approximations 
are the only techniques at disposal, even though one 
should bear in mind their limitations as well as their 
merits. 

Let us consider a general model described by 
F{G,R), satisfying the condition l|H7|) and with 
Ae/y > 0. Thus, we may have a de Sitter Euclidean 
instanton. In the Euclidean version, the associated 
manifold is S4. 

Making use of the instanton approach, for the Eu- 
clidean partition function we have 



Z ~ Z{S4) 

= Z(^\Si)e-^^^^\ (74) 

where / is the classical action and Z*^^) the quantum 
correction, typically a ratio of functional determi- 
nants. The classical action can be easily evaluated 
and reads 



1927r^ 

I{Si) = T^^(Go, Ro) 



(75) 



At this point we make a brief disgression regarding 
the entropy of the above black hole solution. To this 
aim, we follow the arguments reported in Ref. '24!. 
If one make use of the Noether charge method 23 
for evaluating the entropy associated with black hole 
solutions with constant Gauss-Bonnet and Ricci in- 
variants in modified gravity models, a direct compu- 
tation gives 

S — — (^Fro(Go, Rq) + -^Fgo(Go, i?o) 

. ' 

In the above equations. Ah = 47rr|f, rn being the 
radius of the event horizon or cosmological horizon 
related to a black hole solution. This turns out to 
be model dependent. 

Another consequence, as is well known, is the 
modification of the "Area Law" , which reads instead 



2-kA 



H 



A 



H 



4G 



N 



(77) 



One should also stress that, since the above entropy 
formula depends on Fq{G, R), there is always an in- 
determination: any linear term in G appearing in 
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the classical action is irrelevant as far as the equa- 
tions of motion are concerned, while in the entropy 
formula it gives a constant non vanishing contribu- 
tion. This is a kind of indetermination associated 
with the Noether method 25] . 

Furthermore, as in principle the quantity which 
modifies in a non trivial way the usual Area Law 

Fro{Go, Rq) + —FcoiGo, Ro) (78) 

might be negative, there exists the possibilty of 
having negative BH entropies, according to specific 
choices of F{G, R). 

Let us consider an example. For the model defined 
hy F{G,R)=R + .fiG), 



l + f/'(G„)^2(l + :^ 



(79) 



Thus, with the choice /(G) — —aG^, a > 0, one has 
AttAh 1 



S = 



k2 /3- 1 



(80) 



As a result, modulo the Noether charge method in- 
determination, the entropy may be negative. 

Coming back to the general case, we recall that we 
are interested in the de Sitter metric, which reads 



l-^]dt^ 



1 - 



+ r^dSl , (81) 



with Aeff = The Ricci scalar is Rq — 4A. 
Since rn — I, one has 



A 



H 



Yl-K 487r 



i?0 



A, 



(82) 



As a consequence. 



5(54) — f -Fflo(Go, -Rq) + — (Go,i?o, 

(83) 

Taking Eqs. (|67(l and (|75|l into account, from the 
last equation, one obtains 

S{S^) = ^I[S^) , (84) 

which is a good check of our entropy formula (|76|) . 



symmetric space. One-loop contributions are cer- 
tainly important durin g th e inflationary phase, but 
as it has been shown in [30| , they also provide a pow- 
erful method in order to study the stability of the 
solutions. 

We start by recalling some properties of the clas- 
sical model defined by the the choice F{G,R) = 

^(^R + /(G)) , where now the generic function 

/(G) is supposed to satisfy the "on shell" condition 
168|l . Such condition ensures the existence of a con- 
stant GB invariant, maximally symmetric solution 
of the field equations with (^HJ . 

In accordance with the background field method, 
we now consider the small fluctuations of the flelds 
around the de Sitter manifold. Then, for the arbi- 
trary solutions of the fleld equations we set 



(85) 



and perform a Taylor expansion of the action around 
the de Sitter manifold. Up to second order in h^^, 
we get 



S[h] 




Ro + fo 



C2 



, (86) 



where C2 represents the quadratic contribution in 
the fluctuation fleld hf^,y and, in contrast with pre- 
vious sections, here represents the covariant 
derivative in the unperturbed metric goni^- For the 
sake of simplicity, we have also used the notation 

/o = /(Go), A = /'(Go) and h = /"(Go). 

For technical reasons it is convenient to carry out 
the standard expansion of the tensor fleld hfj_^ in 
irreducible components 26], that is 

+\7^\/^a + ^gf,^{h- Aa) , (87) 

where cr is the scalar component, while and 
are the vector and tensor components, with the 
properties 



IV. ONE-LOOP QUANTIZATION OF 
MODIFIED GAUSS-BONNET GRAVITY ON 
DE SITTER SPACE 

Here we discuss the one-loop quantization of the 
class of models we are dealing with, on a maximally 



v^/i(; = 



In terms of the irreducible components of the /i^^ 
fleld, the quadratic part of the Lagrangian density, 
disregarding total derivatives, reads 



C2=Ct 



~\~ ^vector H~ ■^scalar ? 



(89) 
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where Ctensor, elector and Cscaiar represent the 
tensor, vector and scalar contributions respectively. 
They have the form 



4 




Cm 



/o Rq 



h^" , (90) 



Rn 31/ii?g /2i?^ 



'■0 

16 



288 



432 



/o , Ro 17/i-Ro , hRo 



36 



108 



32 



Ai 



(91) 



C 



scalar 



16 48 72 



3 , 5/ii?o , f2Rl 



32 



+ 



32 



hRp 
32 



24 



Ao 



A 



16 



32 



288 




+ - 



3 

32 



/i Rq 
32 

32 



f2Rl 



A2 
^0 



A 



Aqo- 



+h 



Rq fi R^ /2 ^0 
16 



16 



72 



3 3 /i i?o /2 -Rq 



16 



16 



16 



Ao 



hRp 
16 



A 



Aoa, (92) 



where Aq, Ai and A 2 are the Laplace-Beltrami 
operators acting on scalars, transverse vector and 



traceless-transverse tensor fields, respectively. The 
latter expression is valid off-shell, that is, for an ar- 
bitrary choice of the function /(Go). 

As is well-known, invariance under diffeomor- 
phisms renders the operator related to the latter 
quadratic form not invertible in the {h, a) sector. 
One needs a gauge fixing term and a correspond- 
ing ghost compensating term. We can use the same 
class of gauge conditions chosen in Ref. and for 
this reason we refere the reader to that paper, the 
gauge-fixing Cgf and ghost Cgh contributions to the 
quadratic Lagrangian C = C2+ Cgf + Cgh being the 
very same, that is 



C 



af 



a 
2 



e Ai 



i?o 



+ |/^(ao + |^) Aoa 
->i-hAoh--<j ( Ao + ^ ) Aoa 



2 



e [A, + ^] A16 



-f /.(Ao + f) (Ao + f ) Ao. 
-g.(Ao.f)Ao. 



■i6"l^° + T 



Cgh = a 



X ( Ao-f ^ ) Aoo- 



(93) 



+/3 



Rq 



p-3 
AiC^ 



Aqc 



b Ao 



Ro 



X Ao 



Ro 



Aqc 



(94) 



where a,(3,p are arbitrary parameters and C'^,c, 
B^^h are the irreducible components of ghost and 
anti-ghost fields. 

Now, via standard path integral quantization and 
zeta- function regularization ,26- ..27.. .28j . one can 
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compute the one-loop effective action as a determi- 
nant of a differential operator, exactly in the same 
way as it has been done in Ref. Here the tech- 
nical difficulty comes from the fact that the determi- 
nant which gives the one-loop effective action, also 
after simplifications, is a polynomial of fifth order in 
the Laplace operator and thus, in general, it is not 
possible to write it as a product of determinants of 
Laplace-like operators. The general structure of the 
one-loop effective action F^^^ is of the form 



= det(Ao-i?o)det(^Ai + :^ 

-1/2 



X det 



X det[(Ai -f r)(Ai + 

X det(A2 -t-X)-i/2 (95) 

where X, Y, Z, Wk are complicated functions of 
/o,/i,/2- In principle, they can be exactly com- 
puted, but we will not write them explicitly here 
since, as they stand, they are not really convenient 
for direct applications to the dark energy problem. 
However, Eq. (|95|) is certainly useful for a numeri- 
cal analysis of the models, since the eigenvalues of 
Laplace- like operators on de Sitter manifolds are ex- 
actly known and, as a consequence, the one-loop ef- 
fective action can be obtained in closed form. They 
could be also useful for the study of the stability of 
the de Sitter solutions of GB modified gravity, as it 
happens for the simpler case of J-{R) models ^3^]. 

One can see that the structure of Eq. H95|l . here 
written for the gauge parameter /3 = 0, is similar to 
the one for the analog equation (3.24) in Ref. |30l |. 
the only difference being due to the fact that (|95|l 
contains some additional contributions, since the 
starting classical action here depends not only on 
curvature, but also on other invariants. 



V. THE TRANSITION FROM THE 
DECELERATION TO THE ACCELERATION 
ERA IN MODIFIED 
CURVATURE-GAUSS-BONNET GRAVITY 

It is interesting to study late-time cosmology in 
generalized theories, which include both the func- 
tional dependence from curvature as well as from the 
Gauss-Bonnet term (for some investigation of the re- 
lated Ricci-squared gravity cosmology, see [sj). Our 
starting action is ((Tn|l . In the case pm = 0, we can 



consider the situation where F(G, R) has the follow- 
ing form, as an explicitly solvable example: 



F{G,R)=Rf[^ 



(96) 



If we assume that H = ho/t, with a constant ho, 
Eq. (|15|l reduces to an algebraic equation: 




C 



2 (3/io - 1) f'{C) + 3 i2ho - if f{C) , 
2ho {ho - 1) 



(97) 



3(2/10-1)' 

As a further example, we consider the model where 



/ 



Here, it may be shown that 



^ (A, 
3 

'2^ 



2/o 



i?2 



2/o - 2 — - 2/( 



(98) 



(99) 



When /o < and /o > S/Sk^, Eq. ^ has the 
following solutions: 



ho = 



2/o±y^8/o(/o-g^) 



One can solve (I99II with respect to k^/q: 



«'/o = Giho) 



3(2/10-1)' 
'4(/i2 + 2/10-1) 



(100) 



(101) 



and, therefore, by properly choosing /o, we obtain a 
theory for any specified ho- It is easy to check that 
G{oo) = 3, G(l/2) = 0, and G(-1±V2) = oo. 
Since 



G'iho) 



3 {2ho - 1) (3/io - 1) 
2{hl + 2ho-lf 



(102) 



we also get that G{ho) has extrema for ho = 1/2, 
1/3. Moreover, G(l/3) = 3/8. The qualitative be- 
havior of K^/o — G{ho) is given in Fig.^ Then, for 
K^/o < —3, there is a solution describing a phan- 
tom with ho < —1 — \/2 and a solution describ- 
ing effective matter with ho > — 1 -t- \/2. When 
—3 < K^fo < 0, there are two solutions, describ- 
ing effective matter with ho > —1-1- V2. When 
< K^fo < 3/8, there is no solution. When 
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Eq. (|105|l has a real solution for go < 0: 

1/6 



Hn = <- 



12(^2/0 + 3) 



(105) 



Then, the decelerating universe can indeed turn to 
a de Sitter universe, which is accelerating. There- 
fore, the model H103|) could perfectly describe the 
transition from deceleration to acceleration. More 
complicated, dark energy cosmologies may be con- 
structed in frames of such theory. 



VI. HIERARCHY PROBLEM IN MODIFIED 
GRAVITY 



FIG. 1: The qualitative behavior of G{ho). 



3/8 < K^fo < 3/2, there are two solutions describ- 
ing matter with < ho < -I + \/2 < 1. When 

fo > 3/2, there is a solution describing a phantom 
era with —1 — < ho < and a solution describing 
an effective matter era with <l/3</io<— 1 + \/2- 
As one sees, there can be indeed solutions describing 
an effective phantom era, in general. 

Observational data hint towards the fact that the 
deceleration of the universe turned into acceleration 
about five billion years ago. We now investigate if 
we can construct a model describing the transition 
from the deceleration phase to the acceleration one, 
in the present formulation. To this end, we consider 
the following model 



F{G,R) = R 



50 
G 



(103) 



When the curvature is large, as in the case of the 
primordial universe, the last term can be neglected 
and the model I|1U3I) reduces to Then, with 

the choice 3/8 < k^/q < 3/2, it follows < ho < 
— y/2 < 1 and therefore the universe is decel- 
erating. If the curvature becomes small, as in the 
present universe, the last term becomes large. But, 
by including the last term, we have a solution de- 
scribing a de Sitter universe in H17|l . which has the 
form: 



0--5(-Vo + 3) 



90 



12i/4 



(104) 



As /v^/o -I- 3 is positive when 3/8 < n^fo < 3/2, 



Recently the hierarchy problem has been investi- 
gated, in |32j, by using a scalar-tensor theory. Here, 
we will give a somewhat similar but seemingly more 
natural models in the scalar-tensor family with mod- 
ified gravity. In in order to generate the hier- 
archy, a small scale, which is the vacuum decay rate 
Tyac, was considered. Instead of Tyac, we here use 
the age of the universe, ~ 10"^'^ eV, as the small 
mass scale. 

The following scalar-tensor theory [s^ can be con- 
sidered, as an example 

S ^ J d'^x^c"'*'[R~]^d^,(j)d^'(t) 

-Vbe-2'^/'^o~ 



,d^xd^X - U{x) im) 



where a, Vq, and (po are constant parameters and 
U{x) is the potential for x- As the matter scalar, 
X, does not couple with (p directly, the equivalence 
principle is not violated, although the effective grav- 
itational coupling depends on (j) as 

k ^ ne-^ . (107) 
When X = 0, the following FRW solution exists 



a{t) = ao 



= 00 In 



— 2a(4 



2 (a0o + 2) 



{acPo + 2Y 



(108) 



' {a^ + ^) (24«202 , 4a0o + 90^ - 4) 
2Vo 
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The effective EoS parameter is 



Wee 



4 {a^o + 2) 



(109) 



which can be less than —1, in general. 

We may perfectly assmxie the dimensionful param- 
eters K and Vq, and therefore to in H1U6|) . could be the 
scale of the weak interaction ~ 10^ GeV = 10^^ eV. 
If t is of the order of the age of the universe, 
- 10-33 eV, Eq. PUSI) gives 



-a(t>/2 



Then, if 



17 
22 ' 



(110) 



(111) 



K (|107|l is of the order of the Planck length 
(lOi^GeV)" . Therefore, using a model whose ac- 
tion is given by (|106|l . the important hierarchy prob- 
lem might be solved. By substituting Ulllf) into 
(|109|l . one obtains 



61 



33 (0n 24^2) 



(112) 



It is seen that Wqs can be less than — 1 if < 
85/242, and to is given by 



22^ 
612 



\ 



(21^+3) 



2Vo 



(113) 



which is real and positive, as far as Vq > 0. Similar 
method can be applied in the solution of the hier- 
archy problem in a generalized scalar-tensor theory 
including a non-minimal coupling with the curvature 

We now start from the following action for modi- 
fied gravity coupled with matter 



5 = 



(114) 



JF(i?) being some arbitrary function. Introducing 
the auxiliary fields, A and B, one can rewrite the 
action H114|l as follows: 



S 



d X\J—g 



{B iR-A)+ J^{A)} + L 

matter 

(115) 



One is able to eliminate B, and obtain 

-r 1 



S 



d Xy/—g 

~t~ -^matter 



■{T'{A) {R-A)+HA)} 



(116) 



and by using the conformal transformation 

gfi^ e'^g,,^ , (117) 

with 

a = - InT'iA) , (118) 

the action (|116|l is rewritten as the Einstein-frame 
action: 



Se = / d'^Xy/^ 



4 (R-lgP'^dpad^a-Via) 



(119) 



Here, 



V{a) = e"e(e-") -e2-.F(e(e-'^)) 
A J^jA) 

r{A) ~ r{Ay ' 



(120) 



and we solve H118(l with respect to A as A = Q{p)- 
Let us assume that the matter Lagrangian density 
'Cmattor coutaius a Higgs-like scalar field ip 



Under the conformal transformation (|117|l . the mat- 
ter Lagrangian density >Cmattcr is transformed as 



matter 



-Ae2>4 + 



By redefining Lp as 



'^matter acquircs the following form 



1. 



,2„(T 



^ 2 \ 4 , 



(122) 
(123) 

(124) 



where the time derivative of a is neglected. Then, 
the massive parameter /i, which determines the weak 
scale, is effectively transformed as 



(125) 
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In principle, ^ can be of the order of the Planck 
scale 10^^ GeV, but if (fl"^ ^ 10"^^ in the present 
universe, ji could be 10^ GeV, which is the scale 
of the weak interaction. Therefore, there is a quite 
natural possibility that the hierarchy problem can be 
solved by using the above version of modified gravity. 
We may consider the model 



TiR) = R + foR° 



(126) 



with constant /q and a. If a < 1, the second term 
dominates, when the curvature is small. Assuming 
that the EoS parameter, w, of matter is constant 
one gets 



a = aot 
ao = 



2a 



3{l + w) ' 



Po 



(-6/10 + 12/io 



2\ Q-l 



{(1 - 2a) (1 - a) - (2 - a)/io} 



(127) 



and, by using ()127|l . we find the effective to be 
given by 



w 



(128) 



Hence, if w is larger than —1 (as for effective 
quintessence or even for a usual ideal fluid with pos- 
itive w), when a is negative, an effective phantom 
phase occurs where WcS is less than —1. Note that 
this is different from the case of pure modified grav- 
ity. 

By using (|118|) and neglecting the first term in 
(jnU, it follows that 



1 



VfoaR^- 



(129) 



In the present universe, R - (10-33 eV) . Assume 
now that /o could be given by the Planck scale ~ 

1019 GeV = 10^8 eV as /o ~ {lO^^ eY)^^ . Then, 
Eq. ifT^ would yield 



V2 ^ ;loS1("-1) 



(130) 



If we furthermore assume that — 17 = 61 (a — 1), we 
find a = 44/61. In that case, ii w > —1, WcB > 
— 1 and the universe is not phantom like, but (|130|l 
hints towards the possibility that modified gravity 
can solve in fact the hierarchy problem. 



Let us write the action of the scalar-tensor theory 
S = J d^xy/^e°"l'(^R-^d^,(j)d''(j) 



-Voe 



-20/00^1 



This action could be regarded as the Jordan frame 
action. Scalar field (p could be identified with the 
Higgs field in the weak(-electromagnetic) interac- 
tion. Now the ratio of the inverse of the the ef- 
fective gravitational coupling ii — Ke~~ H107|) and 
the Higgs mass fj, is given by 



e 2 



(132) 



S = 



Hence, even if both oi 1/k and n are of the order of 
the weak interaction scale, if e^ ~ 10^^, 1/k could 
be of the order of the Planck scale. 

By rescaling the metric and the Higgs scalar (p as 

gilv e °"''g^,v , ^ e = , (133) 
the Einstein frame action is obtained: 

^ j d^x^{R - ^ (1 + 3a^) d^cjid^'cj) 

-Foe"^'^^^^'^"^"') 

+ (^I—p^-X^^ + --- . (134) 

In the Einstein frame, the gravitational coupling k 
is constant but the effective Higgs mass, jl, defined 

by 



/i = e 2 ^ 



(135) 



can be time-dependent. Hence, the ratio of 1/k and 
p, is given by 



ji Kfl 



(136) 
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which is identical with (|132|l . Then even if both of 
1/k and /i are of the order of the Planck scale, if 
e~^ ~ 10^^, p. could be an order of the weak inter- 
action scale. Therefore the solution of the hierarchy 
problem does not essentially depend on the choice of 
frame. 

Nevertheless, note that the cosmological time vari- 
ables in the two frames could be different due to the 
scale transformation l|133(l . as 



dt dt 



-dt 



(137) 



Therefore, the time-intervals are different in the two 
frames. The units of time and length are now defined 
by electromagnetism. Then, the frame where the 
electromagnetic fields do not couple with the scalar 
field (f) could be physically more preferable. Since the 
electromagnetic interaction is a part of the electro- 
weak interaction, the Jordan frame in (|131() should 
be more preferable from the point of view of the 
solution of the hierarchy problem. 

In the case of /(G')-gravity, whose action is given 
by Hl()|l with (|18|l . it is rather difficult to solve the 
hierarchy problem in the same way, since the factor 
in front of the scalar curvature, which should be the 
inverse of the Newton constant, although it is indeed 
a constant, in the above cases this factor depends on 
time. However, when including a term like g{G)R, 
where g{G) is a proper function of the Gauss-Bonnet 
invariant, the effective Newton constant could be- 
come time-dependent and might indeed help solve 
the hierarchy problem. 

Similar mechanism can work also in -F(G, R)- 
gravity Hl()(l . Introducing the auxiliary fields, A, B, 
G, and D, one can rewrite the action H114|) as fol- 
lows: 



S = / d^xV^ -^{B {R-A)+D{G-G) 



+F{A,G)} + C 
One is able to eliminate B and D, and obtain 

-r 1 



(138) 



:{dAF{A. C) (R-A) 



+dcFiA, C) {G-G)+F{A,G)} 

+^matter • (139) 

If a scalar field is deleted by 

e-'' = dAF{A,G) , (140) 



one can solve (I14U|I with respect to A as ^ 
A{(j,C). Then, we obtain 



S = 



r 1 



-{e-'^ iR-A{a,G)) 



+dcF{A{a, G),G) {G ~ G) + F{A{<j, G), G)} 



(141) 



Varying over G, it follows that G = G, which allows 
to eliminate G 



S = 



xy/-g 



-i{e-- {R-A{a,G)) 



+F{A{a, G),G)}+£ 

n 



(142) 



Performing the scale transformation H117|) . we ob- 
tain the Einstein frame action. If we consider mat- 
ter (Higgs) scalar as in H121I) . we can redefine ip as 
in H123I) . The same scenario as in the case of ^{R)- 
gravity is applied, if e"^/^ ^ 10~^^ in the present uni- 
verse, there is a possibility that the hierarchy prob- 
lem can be solved by working in the above version 
of F{G, R) gravity. Hence, we have proven that the 
hierarchy problem can indeed be solved in modified 
gravity which contains a F{G, R) term. 



VII. DISCUSSION 

To summarize, various types of dark energy cos- 
mologies in modified Gauss-Bonnet gravity — which 
can be viewed as inspired by string considerations — 
have been investigated in this paper. We have 
shown, in particular, that effective quintessence, 
phantom and cosmological constant eras can natu- 
rally emerge in this framework, without the need to 
introduce scalar fields of any sort explicitly. Actu- 
ally, the cosmic acceleration we observe may result 
from the expansion of the universe due to the grow- 
ing of the extra terms in the gravitational action 
when the curvature decreases. In addition, with the 
help of several examples, corresponding to explicit 
choices of the function /(G), we have shown that 
the unification of early-time inflation with late-time 
acceleration in those theories occurs also quite nat- 
urally. Moreover, the framework is attractive in the 
sense that it leads to a reasonable behavior in the 
Solar System limit (no corrections to Newton's law, 
no instabilities, no Brans-Dicke problems appear), 
whatever be the particular choice of /(G). Finally, 
the transition from the deceleration to the accelera- 
tion epoch, or from the non-phantom to the phan- 
tom regime — provided the current universe is in its 
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phantom phase! — may both be natural ingredients 
of our theory, without the necessity to invoke any 
sort of exotic matter (quintessence or phantom) with 
an exphcit negative EoS parameter. 

We have also shown in the paper that modified GB 
gravity has de Sitter or SdS BH solutions, for which 
the corresponding entropy has been calculated. It 
has been explicitly demonstrated that our theory 
can be consistently quantized to one-loop order in 
de Sitter space, in the same way as modified ^{R) 
gravity. 

Dark energy cosmologies in a more complicated 
F{G, R) framework can be constructed in a simi- 
lar fashion too. An attempt to address fundamental 
particle physics issues (as the hierarchy problem), as 
resulting from a modification of gravity, has shown 
that some natural solution may possibly be achieved 
in .7^(i?)-gravity, but probably not in /(G)-gravity 
(albeit the case F{G, R) opens again a new possi- 
bility)). In this respect, it may also be of inter- 
est to study other modified gravities, where a non- 
minimal coupling of the sort Ld-f{G), with La being 
some matter Lagrangian which includes also a ki- 
netic term is introduced. In the specific case of a L^- 
!F{R) non-minimal coupling, such terms may be able 
to explain the current dark energy dominance l36j as 
a gravitationally-assisted one. 

The next step should be to fit the specific astro- 
physical predictions of the above theory with cur- 
rent observational data (for a recent summary and 
comparison of such data from various sources, see 
[ssf). which ought to be modified accordingly, as 
most of them are derived under the (often implicit) 
assumption that standard General Relativity is cor- 
rect. One immediate possibility is to study the per- 
turbation structure in close analogy with what has 
been done for JF(i?)-gravity 0, |33|- This will be 
reported elsewhere. 
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APPENDIX A 

Having in mind the importance of spherically sym- 
metric BH solutions in gravity theories, let us con- 
sider the possibility that F{G, i?)-gravity has the 
Schwarzschild black hole solution: 

ds^ = -e^^dt^ + e-^^dr^ + r^d^l , 

e^-- = 1 - ^ , (Al) 
r 

For simplicity, we consider the vacuum case (T^i/ = 
0) and concentrate on the model (|18() . Then, by 
multiplying gf^^ with 112|) for the action 1181) . we 
find 

= ^R + 2f{G)-2f'iG)R^+Af'{G)R^pR'''' 

-2f{G)R^'"'^R^p,r + Af'{G)R^^R^, 

-2 (V2/'(G)) i? + 4 (V^V^/'(G)) R''^' . (A2) 

In case of the Schwarzschild black hole ljAl|l , one has 

R = Rfiu = , 

G — R^vpaR^ ^ — ■ (-^3) 

Then, Ea. l|A2|) is reduced to be 

= /(G) - G/'(G) , (A4) 

which gives 

/(G) = /oG , (A5) 

with a constant /q. Since G is a total derivative, one 
can drop /(G) in with |(TSJ) for (|X5l) . Hence, 
the Schwarzschild black hole geometry is not a so- 
lution for a non-trivial /(G)-gravity which justifies 
our interest for SdS BH in section HI as for the 
spherically symmetric solution of the above /(G) 
theory. Note that a theory of this sort may contain a 
Schwarzschild solution in higher dimensions, where 
G is not a topological invariant (for a recent exam- 
ple, see K^). It should be also stressed that modified 
gravity of the -F(G, R) form containing more compli- 
cated i?-dependent terms might admit the standard 
Schwarzschild black hole solution. 
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